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We study the superconducting properties of a two-dimensional superconductor in the proximity 
to an electronic topological transition (ETT). In contrast to the 3D case, we find that the 
superconducting gap at T = 0, the critical temperature T c , and the impurity scattering rate are 
characterized by a nonmonotonic behavior, with maxima occurring close to the ETT. We derive 
analytical expressions for the value of such maxima both in the s-wave and in the d-wave case. 
Such expressions are in good qualitative agreement with the phenomenological trend recently 
observed for T™ ax as a function of the hopping ratio t' /t across several cuprate compounds. We 
further analyze the effect of an ETT on the Ginzburg-Landau stiffness r\. Instead of vanishing at 
the ETT, as could be expected, thus giving rise to an increase of the fluctuation effects, in the case 
of momentum-independent electron-electron interaction, we find r\ 7^ 0, as a result of an integration 
over the whole Fermi surface. 



PACS numbers: 74.20.-z, 74.62.-c, 74.40.+k 



I. INTRODUCTION 



The unconventional properties of the normal and su- 
perconducting states of several low-dimensional novel 
electronic materials is a source of continuous interest and 
research. Such materials include the high-T c cuprate 
superconductors (HTSC)lil, as well as some organic su- 
perconductors based on doped BEDT-TTF layers, and 
the ruthenates. In these materials, the interplay between 
their reduced dimensionality and the strength of the ef- 
fective electron-electron interaction is believed to be the 
key for the elusive nature of their normal state, as well 
as for the anisotropic gap characterizing their supercon- 
ducting state. 

A feature common to almost all the material classes 
listed above is a quasi-2D dispersion relation, arising 
from their layered structure and stabilized by the ten- 
dency to confined coherence within layers, due to strong 
correlationso. Indeed, flat bands have been observed in 
nearly all hole- and electron-doped superconductors 
the k phase of BEDT-TTF organic superconductor; 
as well as in the noncuprate layered superconductor 
Sr 2 Ru0 4 (Ref. |). Clear evidence for a 2D Fermi sur- 
face changing topology as a function of doping has been 
recently provided by ARPES measurements in LSCOQ. 
In particular, the role of the proximity to an electronic 
topological transition in establishing the unconventional 
properties especially of the cuprates has been very early 
emphasized (see Ref. ^ for a review). Therefore, in the 



following we will bejnainly concerned with the case of 
the high-T c cupratesa. 

The term electronic topological transition (ETT) has 
been proposed in order to describe the phenomena re- 
lated to a change of the connectivity^ number of the 
components of the Fermi surface (FS)EI. Such a tran- 
sition can be driven by several causes such as isotropic 
pressure, anisotropic deformation, and the introduction 
of isovalent impurities. All these influences can be 
parametrized by their effect on the chemical potential 
fi passing through the critical value £ c , corresponding to 
the transition point. Indeed, such a critical point can 
be well defined only in a pure metal at T-i= where 
a true phase transition of order l\ occural3, according 
to Ehrenfest classification. Typical manifestations of an 
ETT consist in cusp-likepanomalies of physical quantities 
such as the specific heatEH, the density of states (DOS), 
and the conductivity, as well as in the appearance of 
asymmetric singularities of the thermal expansion coeffi- 
cient and thermoelectric power in the dependence of all 
these quantities on z = [i — e c . A non-zero tempera- 
ture or the presence of electron scattering results in the 
smearing of these anomalies and, strictly speaking, in 
washing out the notion of a 2^-order phase transition it- 
self. Moreover, the occurrence of an ETT can be masked 
by an intervening structural transition, as could be in- 
duced by external pressure. The effects of an ETT on 
the properties of metals-and alloys have been thoroughly 
investigated as wellrfc3't3. 

In lower dimensional metallic systems, an ETT is char- 
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acterized by yet stronger anomalies. In particular, the 
DOS of a 2D metal increases logarithmically near an 
ETT, instead of displaying a square-root cusp, as in the 
3D caseEj. Therefore, it has been suggested that an ETT 
may be a clue for the understanding of the anomalous su- 
perconducting state of the high-T c cupratesQ. In partic- 
ular, it is well known that the presence of an ETT in the 
spectrum of a 2D superconductor induces a nonmono- 
tonic dependence of. 1 "*" 3 critical temperature on doping 
or applied pressureOLj, in qualitative agreement with the 
available experimental resultstjO. This has to be con- 
trasted with the 3D case, where an ETT only gpes rise 
to a step-like behavior in the z-dependence of XjiS. 

Moreover, it has been proposed that the proximity 
to an ETT may be the origin of the unconventional 
normal state ,pf. [ tlpjC- | HTSC. In particular, .a, marginal 
Fermi liquidliSEd'ElEa or a non Fermi liquicO behavior 
can be naturally derived for a 2D electron system near 
a Van Hove singularity. More generally, it has been ar- 
gued that the anomalous finite-temperature phenomenol- 
ogy of the cuprates stems from the competition of sev- 
eral broken-symmetry states intervening—near one and 
the same quantum-critical point (QCP)Ejuj. Recently, 
Onufrieva et alEMB showed that an ETT occurring in 
a 2D square lattice with hopping beyond nearest neigh- 
bors is a QCP, with two aspects of criticality: the first 
is related to the singular behavior of the thermodynamic 
properties (Van Hove singularity), while the second is re- 
lated to the existence of the jcritical line T = 0, z > 
of static Kohn singularitiesESO. The proximity to an 
ETT may be characterized by spin density wave (SDW), 
charge density wave (CDW), and d-wave superconduct- 
ing (dSC) instabilities, depending on the appropriate in- 
teraction channels included in the analysis. Onufrieva 
et al. argued that SDW fluctuations dominate in the 
case of the high-T c cuprateso. On the other hand, re- 
cent studies focussed on the competition between AFM 
and AFM-mediated <i-wave pairing via a diagrammatic 
approacb£3, among dSC, AFM, and 7r-triplet pairing 
mean field level in the presence of backward scatteringca, 
among dSC, AFM, and,ppW within the renormalization 
group (RG) approachHElN, among dSC, AFM, and FM 
within the RG and the parquet appr.aachesL3, or between 
dSC and an excitonic ordered stateE3. More recently, it 
has been showntij that elastic umklapp scattering near a 
Van Hove singularity mat, give rise to an RVB-like, in- 
sulating spin liquid stateEj, which exhibits both d-wave 
superconducting and AFM correlations, without being 
characterized by true symmetry breaking, as is typical of 
a quantum ordered state. The competition between su- 
perconductivity and various kinds of density waves in sev- 
eral low-dimensional electron systems in the presence of 
a Van Hove singularity has been reviewed both from the 
experimental and the theoretical point of view in Ref. |3^. 

In this paper, we will concentrate on a single supercon- 
ducting instability (towards either an s-wave or a d-wave 
superconducting state), in the weak coupling limit, thus 
neglecting altogether any other competing ordered phase, 



for a 2D electron system near an ETT. This is of course 
justified only if all other instabilities are characterized 
by weaker couplings, which may not be the case for the 
cuprates. However, such an approximation will enable 
us to derive an analytical expression for the maximum 
gap Ao near the ETT as a function of the band details. 
Such results are in good qualitative agreement with re- 
cent studies of T c in the cuprates correlated with material 
dependent properties, such_as the ratio of next-nearest to 
nearest neighbor hoppingej. 

We will also study the effect of an ETT on the 
Ginzburg-Landau stiffness r\ oc Gi _ , where Gi is the 
Ginzburg-Levanyuk parameter, which characterizes the 
manifestation range of fluctuations near T c (Ref. g). In 
the case of an isotropic FS, r\ is proportional to the square 
of the Fermi velocity. In the vicinity of an ETT, due to 
the presence of 'slow' electrons near the saddle point in 
the electronic spectrum, one may expect an increase of 
fluctuations. However, we will show that, in the case of a 
momentum-independent electron-electron interaction, all 
electronic states on the FS participate in establishing the 
superconducting correlations. Such correlations give rise 
to a superconducting stiffness, whose value is of the same 
order of magnitude p£ the result obtained for an isotropic 
electronic spectrum!^, to the lowest order in z/Ep. 

The paper is organized as follows. In Sec. II, we intro- 
duce a dispersion relation beyond nearest neighbors for 
an electron system in a 2D square lattice, as is typical 
for the cuprates. and discuss its corresponding singular 
DOS. In Sec. [II, we study the superconducting gap A 
at T = as a function of the critical parameter z, in 
the case of s- and (i-wave pairing. In Sec. IV, we dis- 
cuss the effect of impurities on the normal state DOS 
and show that the proximity to an ETT gives rise to a 
nonmonotonic z-dependence of the renormalized quasi- 
particle inverse lifetime t" 1 . In Sec. [v|, we calculate the 
Cooper pair propagator near an ETT, and discuss the 
effects of an ETT on the superconducting fluctuations. 



We eventually summarize in Sec. VI 



II. THE MODEL 



Detailed band structure calc 
density approximation (LDA)t. 



dations within the local- 
j, as well as ARPES ex- 



perimentsEj, show that a realistic tight-binding approxi- 
mation for the band dispersion of most quasi-2D cuprates 
has to be expanded at least up to next-nearest neighbors 
hopping. We then assume the following rigid band dis- 
persion relation for a tetragonal lattice: 

£k = £k — A* = — 2t(cos k x + cos k v ) + it' cos k x cos k y — fj,, 

(1) 

where fi denotes the chemical potential, and the compo- 
nents of the wavevector k are measured in units of the 
inverse lattice spacing. A non-zero value of the hopping 
ratio r — t' /t, measuring the ratio of next-to-nearest vs 
nearest neighbors hopping, slightly modulates the actual 
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shape of the Fermi line £k = 0, and in particular de- 
stroys perfect nesting at fi — as well as electron-hole 
symmetry (Fig. [I]). In order to have a flat minimum in £k 
around the T point, as is. observ ed experimentally for the 
majority of the cupratesE3c3c3, the condition < r < h 
must be fulfilled. The role of an extended saddle point in 
stabilizing superconductivity against other possible low- 
energy instabilities has been established within the renor- 
maliza.tion group (RG) approach in the weak-coupling 
limfO. Moreover, it has been shown that an increase 
of r in the mentioned range correlates with an increase 
of the maximura_T c across different classes of cuprate 
superconductorsE3. However, changes of the shape of 
the FS resulting from screening effects can also corre- 
late with changes in the superconducting properties in 
an indirect way. Indeed, a deformation of the FS also in- 
duces a change of the phase space effectively probed by 
the electron-electron interaction!!]. This is particularly 
relevant for several models proposed for the HTSC, char- 
acterized by effective interactions peaked at X — (0, 7r), 
namely exactly where the Fermi line is most sensible to 
a change in the hopping range r. Possible realizations 
of such a strongly anisotropic k-dependent potential in- 
clude the interaction mediated by antiferromagnetic spin 
fluctuations or by quasicritical stripe fluctuations, due 
to the proximity to a QCP near optimal doping at T = 
(Refs. ^9| , ^0|j5l| ) , as well as electron-electronJjjicractions 
enhanced by interlayer pair-tunneling (ILT)ljH 

As the chemical potential fi in Eq. ([!]) varies from 
the bottom, e±_ = — 4i(l — r), to the top of the band, 
£t = 4i(l + r), the Fermi line £k = evolves from 
an electron-like contour, closed around the T point, to 
a hole-like contour, whose continuation into higher order 
Brillouin zones closes around the M = (7r,7r) point. In 
doing so, an ETT is passed exactly at \x = e c = —it', 
where the Fermi line touches the zone boundaries, and 
assumes the asteroid-like shape depicted in Fig. [I]. It has 
been shown that such a critical value for the Fermi energy 
is stable against tha_RG flow for any repulsive electron- 
electron interactiorO'Ej. Such a result has been recently 
confirmed also when self-energy effects to the quasiparti- 
cle dispersion relation are included, thus demonstrating 
that the pinning of the Fermi surface to a Van Hove sin- 
gularity can actually take place for a rather wide range of 
hole concentration^. The condition fj, = e c corresponds 
to having a saddle point at X = (tt, 0) in the single- 
particle dispersion relation £k, which, for small wavevec- 
tor displacements from X and symmetry related points, 
can be expanded as 




FIG. 1: Fermi line £ k = 0, Eq. ([[]), for a value of the hop- 
ping ratio r = 0.45. As the chemical potential varies from the 
bottom to the top of the band, the Fermi line changes topol- 
ogy, evolving from a closed contour around the T point, to 
a contour, whose continuation in the higher order Brillouin 
zones closes around M = (tt, tt). The change of topology 
(ETT) occurs when the Fermi line touches the zone border, 
i.e. at X — (tt, 0), and symmetry related points. The thicker 
solid lines evidence the hyperbola-like shape of the Fermi line, 
Eq. (0), around X. 



<i£!k|Vk£k| 1 can be computed analytically as 



v(e) 



TT 2 ^fW 



-.K 



£ c £ 



16i 2 - (e + e c 



At 2 



, (3) 



where K (k) denotes the complete elliptic integral of first 
kind of modulus k (Ref. ^). At e — e c , v(e) diverges 
logarithmically. Making use of the appropriate asymp- 
totic expansion for K(k) (Ref. |6C| ), it is then customary 
to identify a regular and a singular contribution to v(e) 
as 



v(e) = v (e) +6v(e), 



(4) 



with vq (e) being a continuous function of e over the whole 
bandwidth such that vq{£ c ) — 0, and 



2 ? 

Pi _ J>2_ = 
2m i 2r7i2 p 



(2) 



where p\ — k x , P2 — k y — tt, and 777,1.2 = [2£(1 ±,2r)] _1 
are the eigenvalues of the effective mass tensortll. Here 
and below, we choose units such that H and the lattice 
spacings are set equal to unity. 

For £_l < £ < £t, the density of states v(e) — 



8v{e) = 2p\n 



4V2/tt 2 



(5) 



where p^ 1 = 47r 2 i\/l — 4r 2 (Fig. ||). In the following, we 
shall often make use of the critical parameter z = /j — e c , 
measuring the distance of the chemical potential from the 
ETT. 
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FIG. 2: Total density of states (v), Eq. ([]), and regular (yo) 
and singular (8v) contributions to the DOS, Eq. (0), as a 
function of energy e, ranging from the bottom (ex) to the 
top (ey) of the band. The inset shows the effect of a non-zero 
energy broadening F (here, T ~ 0.5 % of the total bandwidth) 
on the same quantities [see Eq. (EJ)]. 



III. THE EFFECT OF AN ETT ON THE 
SUPERCONDUCTING GAP 

We will now discuss the effects of the proximity to an 
ETT on the superconducting properties of a 2D, single- 
layer system, both for an s-wave and for a d-wave order 
parameter, in the weak coupling limit. Our starting point 
will be the BCS equation for the gap function Ak, which 
at T = reads as 

Ak = --Vfe — ■ (6) 

N v 2E *' 



Here, = + | Ak| 2 is the upper branch of the su- 
perconducting excitation spectrum, Vkk' denotes the in- 
terparticle potential, and the sum runs over all the N k 
points in the 1BZ. In the case of s-wave symmetry, we 
assume Vkk' = —X, i.e. a constant over the whole 1BZ, 
whereas in the cZ-wave case we take the potential in the 
separable form Vkk' = —Xgugw, ffk = ^(cosk x — cosk y ) 
being the lowest-order lattice harmonic corresponding to 
d-wave symmetry. Accordingly, one has Ak = Ao in the 
s-wave case, and Ak = Ao<?k in the d-wave case, respec- 
tively. It is worth emphasizing that, in both cases, the 
coupling constant A > has been assumed independent 
of doping. This amounts to neglecting higher-order cor- 
relation effects among interacting particleaia. Moreover, 
the weak coupling hypothesis allows us to neglect renor- 
malizations of the shape of the Fermi surface, which are 
certainly expected in the strong coupling limit, and are 
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known to give rise to another kind of ETT as welM 

Eq. (||) implicitly neglects the possibility of any pair- 
ing instability other than singlet superconductivity in the 
Cooper channel (characterized by a pair relative momen- 
tum P = 0). Possible alternative intervening pairing 
instabilities include- e. 17., antiferromagnetism and the tt- 
triplet paired stateEl Such instabilities would be charac- 
terized by large momentum transfer near the 'hot spots' 
(0, 71") and (7T,0). They have been shown to coexist and 
win out singlet superconductivity at a mean field level 
near hatt-filling, when backward scattering is a relevant 
processE2l. However, within our weak-coupling approx- 
imation, it is consistent to retain only one kind of in- 
stability (namely, Cooper pairing in the singlet, P = 
channel, with either s- or (i-wave symmetry), under the 
assumption that other instabilities are characterized by 
weaker couplings. 

We will first analyze the gap equation, Eq. (||), close 
to an ETT (\z\ <§; 4i). In the s-wave case, the sum- 
mation over the 1BZ in Eq. (^) can be transformed into 
an integral over energy weighted by the DOS, which we 
approximate by its singular part 5v[e) in Eq. (^). The 
Fermi line corresponding to the ETT divides the 1BZ in 
two regions, £k < and £k > 0, which are electron-hole 
conjugated of each other. Separating the contributions 
coming from these regions, the gap equation can be com- 
pactly written as: 

-L =5+ + 5 _ ; (7) 
Xp 

where S± represent the pairing susceptibility integrated 
between the ETT and either band edges (see App. |A| for 
details). 

While in 3D BCS theory S± are logarithmically diver- 
gent in the limit Ao — * (Ref. [32]), the proximity to an 
ETT in 2D makes them divergent as ~ In 2 Ao. Direct in- 
spection of Eq. (^) as well as numerical calculations show 
that Aq(z) is maximum near the ETT. Such a nonmono- 
tonic dependence of the superconducting gap on the criti- 
cal parameter z is in agreement with the phenomenology 
of the HTSC, where Ao and T c arep-eharacterized by a 
parabola-like dependence on dopingta. This has to be 
contrasted with the step-like behavior observed in the 
3D casetLS, where it has to be emphasized that no diver- 
gence occurs in the DOS at the ETT. However, due to 
the electron-hole symmetry breaking induced by a non- 
zero hopping ratio r, Aq(z) is not an even function of 
z, and its maximum will actually occur not exactly at 
z = 0, as will be discussed below. Solving Eq. (0) for A 
in the weak-coupling limit (Xp <C 1, Ao <C 4t), at z = 
one finds: 



» , s Sv^ / H 1, 2 l + 2r \ , 

A (z = 0) * -5- exp + 5 In 2 — J , (s-wave) (8) 
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which can be then taken as a first approximation to the gap maximum at T = 0, in the s-wave case. Following the 
same procedure, qualitatively similar results can be derived for the critical temperature T c as a function of z (see also 
Refs - & 

In the d-wave case, due to the anisotropic k-dependence of the integrand in Eq. (|6|), it is not possible to explicitly 
separate the integration over energy, and a different approach must be followed (see App. ^ for details). However, 
the proximity to an ETT does endow the pairing susceptibility with an analogous asymptotic low-Ao behavior, like 
in the s-wave case, which eventually results in the following weak-coupling expression for the gap amplitude at T = 0, 
z = (A/3< 1, A < At): 



A„(: =0)^ -l//iir! o.xp( - v /-L + / 2 ( r ) ) . !,Mvav<0 (!)) 



wher e /i(r ) = 2&~ 1 (1 - 4r 2 )(Vl + 2r + y/l - 2r)- x , f 2 (r) = ln 2 (&7T\/l - 4r 2 ) + 
2 In v 7 ! - 4r 2 ln[27r _ Vl - 4r 2 (vT + 2r + vT^r)" 1 ], and 6 = e 2 /8. 



Fig. || shows Ao(0) both in the s- and in the d-wave 
case, Eq. and (Q), respectively, as a function of the 
hopping ratio r, for several values of A/t. In view of the 
fact that T c oc Ao, as in any mean- field theory, Fig . is in 



39, show- 



good qualitative agreement with Fig. 5 of Ref. 
ing a direct correlation between the experimental T™ ax 
and the hopping range r for several cuprate compounds. 
Moreover, our results suggest that such an effect is a 
general consequence of the proximity to an ETT, and 
is roughly independent of the superconducting pairing 
symmetry. 
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FIG. 3: Normalized gap amplitude Aq(z — 0,r)/Ao(z — 
0, r — 0) at T — 0, as a function of the hopping ratio r = t'/t, 
for different couplings A/t = 0.9 -f- 1.1 (bottom to top). Con- 
tinuous lines refer to the s-wave case, Eq. (|^), while dashed 
lines refer to the d-wave case, Eq. (bl). One can recognize 
the direct correlation between T"™ ax oc Ao(z = 0) and r, as 
observed in Ref. |39| . 

Expanding Aq(z), as implicitly defined by Eq. (Q) 
around z — 0, one finds that the maximum of Ao ac- 
tually occurs at a larger value of the critical parameter, 
which, in the s-wave case, is given by 



(10) 



A qualitatively analogous result applies to the d-wave 
case. Therefore, as an effect of the electron-hole asym- 



metry induced by a non-zero hopping ratio r, the max- 
imum in Ao is actually located in the hole-like region 
(z max > 0), in agreement with the phenomenology of 
some hole-doped cuprate compounds. For instance, a 
representative high-T c cuprate such as LSCO is charac- 
terized by an optimal doping level of x op t — 0.15, ly- 
ing in the hole-doped region, while a doping-dependent 
crossover from a hole-like to an electron-like FS has been 
clearly observed at a somewhat larger doping x c ~ 0.20 
(Ref. ||). On the contrary, no direct evidence of a change 
in the FS topology has been so far reported for Bi-2212 
(Ref. |6^), whose FS displays a hole- like character at all 
dopings, including optimal doping. This implies that the 
ETT is located at a much larger distance from optimal 
doping, which is consistent with Eq. ([To]) above, given the 
larger value of the gap amplitude of Bi-2212 than that of 
LSCO. 



IV. EFFECT OF IMPURITIES 

We now turn to consider the more realistic case, in 
which electron scattering from non-magnetic impurities 
is included. Here, we will be mainly concerned with 
the normal state properties. A finite quasiparticle life- 
time induces a broadening of the energy linewidth of a 
quasiparticle state. Therefore, the use of a quasiparti- 
cle description and the definition of a Fermi surface for 
impure metals can, at first sight, be objected. Indeed, 
quasimomentum is a 'good' quantum number only for 
electrons moving in a periodic potentials. Scattering of 
electrons on impurities results in momentum relaxation 
and in the corresponding smearing of the Fermi surface 
in momentum space. The characteristic scale of such a 
smearing is ~ t _1 , where r is the elastic relaxation time 
at low temperatures. The value of r can easily exceed 
the quasiparticle energy ~ T even for moderate impurity 
concentrations. 

Nevertheless, elastic scattering does not result in en- 
ergy relaxation. This means that, in principle, one can 
solve exactly the eigenvalue problem for the Hamiltonian 
of the electron in a lattice with some specific realiza- 
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tion of the impurity potential. The eigenstates of such 
Hamiltonian can then be chosen as a basis in the Hilbcrt 
space and the 'Fermi surface' in this space can be defined 
as the surface separating the low-energy occupied eigen- 
states from the high-energy empty eigenstates at zero 
temperature. It is evident that the Fermi surface defined 
in this way does exist, and that the elastic scattering has 
no effect on the quasiparticle lifetime in the vicinity of 
the Fermi surface (see also Refs. p|.|l2|). 

The effect of a nonvanishing impurity scattering rate 
can then be accounted for in the DOS by means of a 
convolution between v(e) and a Lorentzian of finite width 

r : 



vr{e) 



r 



7T (£ - ef + r 2 



"(0- 



(ii) 



Such a procedure£3 effectively smears out the logarithmic 
singularity in the DOS at the ETT into a pronounced 
maximum of finite width ~ T (see inset in Fig. ||). 
Nonetheless, it is still possible to separate a 'regular' and 
a 'singular' contribution to ^r(e) as 



vr(e) = i4(e)+6vr(e), 
with Svr{s) now given byEl: 



Svr(e) = 2 pin 



V(e-e c ) 2 +r 2 



(12) 



(13) 



From a physical point of view, the energy linewidth 
broadening associated to impurity scattering has the ef- 
fect of 'blurring' the Fermi line. As a consequence, one 
expects that the ETT occurs slightly farther from fi = e c 
(i.e., for \z\ > 0), as soon as such a blurred Fermi line 
touches the border of the 1BZ0. 

Thus far, we have assumed a constant energy linewidth 
r over the whole band. This is clearly an approxima- 
tion, since the quasiparticle lifetime rt is generally an 
anisotropic quantity over the 1BZE3. The last statement 
holds true even in the simplest case of isotropic impurity 
scattering, due to the anisotropy of the single-particle 
band structure. In particular, the proximity to an ETT 
in 2D endows the quasiparticle lifetime with a nonmono- 
tonic behavior, in contrast to, the 3D case, where a step- 
like z-dependence was foundQ. Following Ref. ^, one can 
write the self-consistent equation for the retarded quasi- 
particle self-energy X R due to impurity scattering as 



E R (w,z) 



w J p [ep 



+ Z + LO 



£ R (c,z)] \ 
(14) 



where e p is the asymptotic single-particle dispersion re- 
lation near the saddle point defining the ETT, Eq. (@), 
rij denotes the concentration of impurities, and Uq is the 
impurity scattering strength, here assumed independent 
of p. Performing the integrations as in Refs. plEl, but 



now for the 2D case, one arrives at the self-consistent 
expression: 



£ = — In 

2r 



y/T+Z + v 7 ! + w + £ H 



+ £ R 



(15) 



where t^ 1 = TT~ 1 rii\uo\ 2 (mim2) 1 / 2 , and all energies 
are in units of a cut-off energy ~ p _1 . Fig. ^ shows 
the renormalized quasiparticle inverse lifetime r _1 = 
— 2Im£ R (u; = 0, z) as a function of the critical parameter 
z, for several values of t^ 1 . As anticipated, one observes 
a maximum in r _1 at z > 0, as an effect of crossing an 
ETT. 
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FIG. 4: Renormalized quasiparticle inverse lifetime r _1 = 
— 2ImE H , Eq. (|lE|), resulting from isotropic impurity scatter- 
ing. Different curves correspond to increasing values of Tq 1 
(bottom to top). The proximity to an ETT induces a non- 
monotonic dependence of r _1 on the critical parameter z, 
with t -1 assuming its maximum value at z > 0. 



The study of a 2D superconductor in the dirty limit 
goes beyond the reach of the present analysis. The depen- 
dence of Ao as well as of the DOS on the impurity concen- 
tration has been derived in the d-wave case in Ref. for 
an isotropic dispersion relation. In d-wave superconduc- 
tors, Dirac-like single-particle excitations can be , 
at virtually no energy cost near the gap nodestJOl 
Within the QCP scenario, long-range interaction between 
such gapless modes is mediated via the fluctuations of an 
intervening order parameter at T = 0. Current proposals 
for the HTSC include the possibility of the proximity to 
a quantum ordered phase characterized by either charge 
or (AFM) spin fluctuations, as well as fluctuations re- 
lated to the opening of another subdominant contribu- 
tion to the superconducting OP, usually accompanied by 
time-reversal breakingcJS. Recent results for 2D d-wave 
superconductors in the presence of disorder yield correc- 
tions to the density of states coming bpth from the diffu- 
sion (Q = 0) and the Cooperon modeO, as well as from 
the diffusive mode with Q = (tt, n)E3. 
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V. GINZBURG-LANDAU STIFFNESS NEAR 
AN ETT 

As is well known, the normal state of HTSC is charac- 
terized by several anomalous properties at the transition 
edge. Such properties include a peak in the c-axis re- 
sistivity, an anomalously large sign-changing c-axis mag- 
netoresistance, as well as the opening of a pseudogap, 
which is observed both in the c-axis optical conductivity, 
in tunneling experiments, and in the NMR relaxation 
rate (see Ref. ^ for a review) . On the basis of the Fermi 
liquid theory, it has been recently demonstrated that the 
renormalization of the one-electron DOS in the vicinity 
of the Fermi level due to the electron-electron interac- 
tion in the Cooper channel is able to explain satisfacto- 
rily many of these pseudogap-like manifestations bothjii 
the overdoped and in the optimally doped compoundscJ. 
Moving across the phase diagram of the HTSC from 
the overdoped, bad metallic region, towards underdop- 
ing, the enhancement of the mentioned effects correlates 
with an increase of the Ginzburg-Levanyuk parameter, 
Gi(2D) &iT c /E-p, thus making perturbation theory less 
reliablec2l. Nevertheless, such a rapid growth of the nor- 
mal state anomalies with the decrease of doping strongly 
overcomes the theoretical predictions, thus making it dif- 
ficult to attribute such an effect to the mere shrinking of 
the FS. 

Indeed, ARPES studies indicate a marked increase 
of the FS anisotropy in the a6-plane with underdoping, 
which is accompanied by the development of an extended 
saddle point in the electronic spectrumc3. One can iden- 
tify two characteristic energy scales related with such a 
FS, namely the size of its 'bulk' part, E-p w 0.3 eV, and 
the 'width' of the saddle point, \z\ = |/z — e c \ ~ 0.01 eV. 
The large difference in the magnitude of such energy 
scales is therefore suggestive of a crossover, related to 
the special role played by the electronic states lying close 
to the saddle point (the so-called 'slow' electrons). Intu- 

I 



with Vl v = 2itTv the bosonic Matsubara frequencies. 

The superconducting critical temperature can be char- 
acterized as the temperature at which L presents a pole 
at q = and f2„ = 0. The procedure to deal with the 
integral in Eq. ( p"8| ) is outlined in App. |B| One eventually 
arrives at the result: 

11(0, 0; z,T)=™T £ -lnf-^f-,), (19) 

w n >0 \ n < / 



itively, one may expect a replacement of E-p with a small 
z in the denominator of Gi( 2 D)j which would result in 
the breakdown of the perturbative approach developed in 
Ref. [l(j| Nevertheless, as is demonstrated below, the sole 
existence of an ETT in the electronic spectrum is not able 
to change the character of the isotropic electron-electron 
interaction in the Cooper channel. This implies that the 
breakdown of fluctuation theory in the underdoped com- 
pounds has to be related to some other properties of the 
HTSC. In order to substantiate the above statements, in 
the present section we shall derive the momentum depen- 
dence of the two-particle Green function in the Cooper 
channel near T c and close to an ETT. 

In the case of a 2D electron system characterized by 
the approximate spectrum given by Eq. (|^), close to an 
ETT, the single-particle Green function can be written 
as 

G(p,uj n ;z) = (iu> n - e p + zy 1 . (16) 

Here, the electron quasi- momentum p is measured from 
the saddle point location, as in Eq. (||), and uj n — 
27rT(n+i) are the fermionic Matsubara frequencies. As 
already emphasized in Sec. II, the condition z > de- 
scribes the case of an open Fermi surface without any 
disrupted neck, whereas the opposite one, z < 0, is ap- 
propriate to a closed Fermi surface with respect to the 
r point (Fig. |l]). The two-particle Green function in the 
Cooper channel L(q, Q,^) can be expressed within the lad- 
der approximation by means of the polarization operator 

L-\q, Q v ;z, T) = A" 1 - n(q, Sl v ;z, T). (17) 

Here, A > denotes the momentum independent effec- 
tive electron-electron interaction, q is the Cooper pair 
momentum, and 



(18) 

I 

the saddle point, and the Debye frequency uip> has been 
introduced as a cut-off in the summation over the Mat- 
subara frequencies. The equation for the critical temper- 
ature then reads 

A- = ™T C £ -mf^). (20) 



U(q,n u ;z,T)=Tj2 J T^jTg ^(p + q, w n +„; z)G(—p, —uj n ; z) = T ^ 7(q, ( 



where m = ^Jm\mi is the geometric average mass around 



For small z (\z\ <ti ^d), one recovers the well-known 
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resul 



T 



WD 

2^ 



exp 



(21) 



where p — m/ (2-ri 2 ) denotes the DOS at the saddle point. 
Eq. (El]) is in agreement with the s-wave result for Ao, 
Eq. (M), when the assumption of a phonon- mediated pair- 
ing mechanism is made and the limit r — > is taken. 
Vertex and cross corrections to Eq. ( |2l| ) in terms of the 
Migdal adiabaticity parameter ujd /E-p have been shown 
to decrease_th£ enhancement of T c due to the proximity 
to an ETTLJO. Analogously, for the two-particle Green 
function close to the superconducting transition and in 



J 



the proximity of an ETT (\z\ «T~ T c ), one finds 



L-\0,0;z,T) = -pin 



( _wd_\ T — T c 
\2«T C 



(22) 



One observes that the presence of the ETT results in the 
appearance of the additional large factor ~ ln(o;D /T c ) in 
front of the reduced temperature. 

In order to determine the superconducting stiffness 
tensor rji(z), one is led to consider the q-dependcncc 
of the polarization operator, Eq. (p"8"|). Expanding 
J(q, u> n , — oj n ; z) in Eq. ( |l8| ) for small q and z, up to 
quadratic order in q, one has 



I(<1, u n , z) = I Q {0,w n , -w„; z) + h(q 2 , u„, -u n ; 0) + I 2 (q 2 ,uJn, -w„; z), 

I 



(23) 



where Ij are defined in App. [El, and Iq has been used 
above for the definition of the critical temperature, 
Eq. (|l^). One eventually finds for the q 2 -dependence of 
the two-particle Green function in the limit \z\ <C T ~ T c 
the relatively classical form: 



L-\q,0;0,T) 



In 



2ttT c 



T-T c 
T c 



mil + mil 



(24) 

where the components of the superfluid stiffness tensor 
are given by 

7C(3)£ F 



8n 2 T 2 7 



(25) 



to the lowest order in z/ E-p. The latter expression is anal- 
ogous to the result obtained in the standard 2D isotropic 
case, with £ p = p 2 /(2m) — /i, where the superfluid stiff- 
ness reads r\ = 7((3)E F /(16TT 2 T 2 m), and the DOS per 
spin is Pjzp) = m/(2ir) (Ref. ^J). It is worth noting that 
in Eq. ( |24| ) the effective mass of the fluctuating Cooper 
pair gets increased by a factor ln(uip, / 2ttT c ) , with respect 
to the case in which a parabolic spectrum is assumed. 
This implies a reduced role of fluctuations near an ETT. 
Indeed, our results demonstrate that the temperature 
range of the fluctuation regime is governed by essentially 
the same Gi, while the propagator's effective mass is en- 
hanced. 

Summarizing, the results obtained above shows that a 
topological singularity in the electronic spectrum prac- 
tically does not affect the Ginzburg-Landau_stiffhess, in 
contrast to what was intuitively speculatecO. The rea- 
son thereof is that the value of r\ is formed by gathering 
the contributions of the electronic states belonging to the 
whole Fermi surface, not only by the 'slow' ones. Finally, 
we note that in the approach we followed, only the po- 
larization loop Eq. (n8h is critical, since Cooper pairing 



of non-zero center-of-mass momentum is not taken into 
accounts 



VI. CONCLUSIONS 

We have reviewed the effects of an electronic topologi- 
cal transition on the superconducting properties of a 2D 
electron system, with an energy spectrum characterized 
by a minimum at the T point and an extended, dop- 
ing dependent saddle point at (tt, 0), as is typical for 
most single-layered, hole-doped HTSC. We analytically 
derived the expressions for the superconducting gap Ao 
at T — close to an ETT, both in the s-wave and in 
the d-wave case. In contrast to the 3D resulttl, Ao turns 
out to be characterized by a nonmonotonic behavior as 
a function of the critical parameter z, with a maximum 
at z ~ 0, i.e. close to the ETT. Due to the electron-hole 
symmetry-breaking induced by a non-zero value of the 
hopping ratio r, we find that the maximum of Ao actually 
occurs at z > 0, i.e. it is slightly shifted towards the hole- 
like region, as observed in-XSGQEl. We point out that in 
previous calculationsEjtj'CJ'Lirca the maximum of T c as a 
function of the critical parameter z occurs at the ETT, 
z = 0. This result has been often used as an objection 
against the relevance of the Van Hove scenario for the 
cuprates, since ARPES data show that the optimal dop- 
ing does not correspond to the critical point z = 0, and 
that the FS preserves a hole-like character over the entire 
doping range for almost all hole-doped compounds (see, 
however, also the recent results of Ino et al. for the LSCO 
compounds). On the contrary, we have shown that the 
observed difference between optimal doping and the dop- 
ing actually corresponding to the ETT can be justified 
by taking into account an electronic spectrum beyond 
the hyperbolic approximation. Moreover, we find that 
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the dependence of Ao on the hopping ratio r is in good 
qualitative agreement with the phcnomcnological results 
collected by Pavarini et al. for several HTSC materialsE3, 
thus demonstrating the role of the band structure pecu- 
liarities and, in particular, of next-nearest neighbor hop- 
ping, in stabilizing high-temperature superconductivity 
in the cuprates. 

In the presence of impurities, the Fermi line is effec- 
tively smeared, and one expects the anomalies due to 
the proximity to an ETT to occur at a larger value of 
the critical parameter, as soon as such a 'blurred' Fermi 
line touches the zone border. We also derived the energy 
dependence of the retarded quasiparticle self-energy S R 
due to impurity scattering in the 2D case, for a simpli- 
fied, hyperbolic dispersion relation. In contrast to the 
3D case, S R is again characterized by a nonmonotonic 
^-dependence, thus confirming that the quasiparticle life- 
time Tk is generally an anisotropic quantity over the 1BZ. 

Finally, we addressed the issue of the range of fluctu- 
ations near T c . By explicitly computing the two-particle 
propagator in the Cooper channel near T c , we derived 
the expression for the supcrfiuid stiffness rj close to an 
ETT. Although the Fermi velocity vanishes at the saddle 
point, we find a non-zero value of 77, in complete analogy 
with the Ginzburg-Landau result for an isotropic elec- 
tronic spectrum, thus showing that all electronic states 
participate in establishing the superconducting correla- 
tions. Moreover, our results show that the role of fluc- 
tuations even diminishes near the ETT. Indeed, while 
their temperature range is determined by about the same 
value of Gi, the effective mass of the fluctuating Cooper 
pair increases by a factor ln(a>D/27rT c ) with respect to 
the case of a parabolic spectrum. Therefore, in the de- 
nominator of any fluctuation contribution there appears 
a large logarithm, which implies a relative suppression of 
fluctuation effects near the ETT. 



Nicosia for warm hospitality during the period in which 
the present work was brought to completion. 



APPENDIX A: EVALUATION OF THE PAIRING 
SUSCEPTIBILITY IN THE d-WAVE CASE 

We briefly outline the derivation of the asymptotic de- 
pendence of the integrated pairing susceptibility close to 
an ETT in 2D. In the s-wave case, changing integration 
variables in the gap equation, Eq. (^|), from wavevector 
k to energy, we can write the pairing susceptibility inte- 
grated between the ETT and either band edges as 



S'4 



ln(«±£K 

" -/(£tC±) 2 + 4 



(Al) 



where we have introduced the dimensionless auxiliary 
quantities C± = z/{At ± 8t'), 5± = A /(4i ± 8i'), and 
a± = (1 ± 2r)/(4y2\/l-4r 2 ). 

In the d-wave case, the double integration over 
wavevector k cannot be reduced to a simple integra- 
tion over energy, and one has to change variables to 
£ = ±(e k -e c )/[4t(l±2r)] > g = gu- In such case, Eq. (§) 
can be written as 



~2A~ 



D_i_ + D_ , 



(A2) 



where the integrated pairing susceptibility now reads: 



D± = 4 



dg 



(A3) 



Acknowledgments 

We thank J. V. Alvarez, G. Balestrino, B. Ginatempo, 

F. Onufrieva, P. Pfeuty, P. Podio-Guidugli, R. Pucci, and 

G. Savona for useful discussions. Partial support from 
MURST through COFIN fund 2001 and the E.U. through 
the F.S.E. Program (G. G. N. A. and E. P.), and from 
I.N.F.M. and I.S.I. (Torino) through the 'Progetto Stage' 
(G. G. N. A.) is gratefully acknowledged. G. G. N. A. also 
acknowledges the D.S.T.F.E. "Tor Vergata" and P. G. 



and 



Ji = (l-2rg) 2 +4r[g-r + (l±2r)£], (A4a) 
J 2 = (l+ 5 ) 2 -(l- v/ji) 2 /(4r 2 ), (A4b) 
J 3 = (1 - gf - (1 - v^i)7(4r 2 ). (A4c) 



Eq. (A3) leads to hyperelliptic integrals, that cannot Jae 
generally expressed in terms of known special functionsLj. 
We next set il = -(1 - \/Jx)/{2r), with ft -> £ as r -> 0, 
with which D + in Eq. (|A3j) transforms into 



- I dn 



dg 



g 



1 



^ [r{v 2 _ g 2) _ z] 2 + A 2 g 2 + n) 2 
I 



g 2 ][{\ - n) 2 - g 2 



(A5) 



where all energies are in units of At and the integration is be to performed over the curvilinear triangle defined by 
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9>0,g< 1-0, and v 2 -g 2 > 0, with v 2 = l + n 2 + fl/r. 
Such triangle represents the hole-like region of the 1BZ, 
£k > £o hi these new coordinates. 

In the limit r — > 0, one branch of the hyperbola defined 
by = t> reduces to the fi = axis, and r(v 2 — g 2 ) — ► Q. 
A further change to polar coordinates as O = pcos9, 
g = psm9 then a llow s to exactly decouple the two in- 
tegrations in Eq. (|A5|), the integration over p leading to 
elliptic integrals. Extracting the logarithmic divergence 
of these latter near 9 — tt/2 (corresponding to the ETT, 



together with p = 1) yields the answer [see Eq. ( |Aq ) 
below, with r = 0]. 

In the case r ^ 0, we could not find any such simple 
change of variable s, a llowing to exactly decouple the in- 
tegrations in Eq. (A5). However, since only the behavior 
of the dispersion relation close to the ETT is believed 
to determine the asymptotic properties of the integrated 
pairing susceptibility, we may linearly expand r(v 2 — g 2 ) 
near Q, = 0, g = 1 and set f2+2r(g — 1) = pcos9. Within 
such approximation, Eq. (|A5|) then reads: 



7r/2 



dO 



\a-P+\ 



cos z 



An sin 



2 Q 1 - 4r 2 



dp 



+ p)(a- - p){f3+ + p)(0- - p) ' 



(A6) 



where 



l±2r 



cos6> 



- (1 - 2r) sin 9' 
l±2r 



cos (9- (l + 2r) sin 9 



(A7a) 
(A7b) 



Following standard methods (see, e.g., Ref. |7^), the inner 
integral can be now expressed as a combination of ellip- 
tic integrals, which diverge logarithmically as 9 — > tt/2, 
whence Eq. (^) follows. 



APPENDIX B: EVALUATION OF THE 
POLARIZATION OPERATOR 



The momentum integration in Eq. ( |l8|) for the polar- 
ization operator can be performed by reducing the inte- 
gration domain to the first quarter of the 1BZ, and divid- 
ing the latter into the two triangles defined by {p : p 2 < 
(m/mi)pi; p\ > 0}, and {p : p 2 > {m/mi)p\; pi > 0}, 
respectively. One finds: 



II(0,0;z,T) = T^ri(p,w n ,—u n ;z) = ^ [/(w,,, -w n ; z) + /(w„, -w n ; -z)] 



(Bl) 



where 



>o Jo 

Performing the inner integration, one obtains 



f(ui n) -uj n ;z)= / dxi I dx 



„2 „2 



2 „2 



X 2 ~ X l ~ lUJ n + Z XT, — Xi + VjJ n + Z 



2 ^1 



f(u ni -uj n ;z) 



1 



2«w„ Jo 



dx 



ln(x — iV z — iuj n — x 2 ) ln(a; + i\[z — iu n — x 2 ) 



Z — lOJ n — X* 



H.c. 



Z — llVn — X* 



(B2) 



(B3) 



The last integration can be carried out observing that we obtain: 



ln(x ± i\Jz - iuj n - x 2 ) t d 2 / : ~ 

, . =z = ±773- In {x±i\/ Z- lUJ n - ar), 

\J z — iu) n — x z <i ax 



Tin \ m 1 1 W D 
1(0, u n , -w„; z) = — rln- 



, , , - , ,• (B5) 
2tt \u) n \ z z +u>^ 

In order to calculate the GL stiffness, let us now ex- 
(B4) pand the polarization operator up to quadratic order in 
q. One has: 



and choosing the branches of the logarithms in order to 
make them complex conjugated of each other. Finally, 



n(g, 0; z, T) = 11(0, 0; z, T) + p £ m 2 , (B6) 
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where the components r\i of the stiffness tensor are de- 
fined by 



P E WVi = T^ihiq 2 ,uj n , -uj n ;0) + I 2 (q 2 ,iv n ,-u) n ; z) + . . .]. 

i=l,2 oj„ 

One explicitly finds: 

Vl{z) = 2^'^j x\-x\- iw n • 1 - - 



(B7) 



{x\ — x\ + iuj n ) 3 2m 



9m, 



{x\ -x\ + iuj n ) 3+k {x\ - x\ - iLu n ) 2 - k 



(B8) 



to the first order in z/Ep. The first integral, giving the principal contribution to the stiffness, can be evaluated using 
polar coordinates as: 



2ir r co „3/„4„„2 



rp ^ pZTT POO 

77,(2 = 0) = _ V — / dxj> I dr 
mi f-^ uj 2 Jo Jo 



r 3 (r 4 cos 2 2(f> — 1) cos <j> 
(r 4 cos 2 20+ l) 3 



ttT x ^ 1 

mi^ W 3 F 8 7r 2 TW 



/>oo 

/ dt 

la 


' d 


d 2 - 


1 




da 


da 2 _ 


yfOLyJt + a 


a=l 



7C(3)^f 



(B9) 
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